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ABSTRACT 

We present an accurate analytic fitting formula for the numerical results for the relativistic 
corrections to the thermal Sunyaev-Zel'dovich effect for clusters of galaxies. The numerical 
results for the relativistic corrections have been obtained by numerical integration of the collision 
term of the Boltzmann equation. The fitting is carried out for the ranges 0.02 < 9 e < 0.05 
and < X < 20, where e = kBT e /m e c 2 , X = Huj/kBTo, T e is the electron temperature, to 
is the angular frequency of the photon, and To is the temperature of the cosmic microwave 
background radiation. The accuracy of the fitting is generally better than 0.1%. The present 
analytic fitting formula will be useful for the analyses of the thermal Sunyaev-Zel'dovich effect 
for high-temperature galaxy clusters. 

Subject headings: cosmic microwave background — cosmology: theory — galaxies: clusters: 
general — radiation mechanisms: thermal — relativity 
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1. INTRODUCTION 

Compton scattering of the cosmic microwave background (CMB) radiation by hot intracluster gas 
- the Sunyaev-Zel'dovich effect (Zel'dovich & Sunyaev 1969; Sunyaev & Zel'dovich 1972, 1980a, 1980b, 
1981) — provides a useful method to measure the Hubble constant H (Gunn 1978; Silk & White 1978; 
Birkinshaw 1979; Cavaliere, Danese, & De Zotti 1979; Birkinshaw, Hughes, & Arnaud 1991; Birkinshaw 
& Hughes 1994; Myers et al. 1995; Herbig et al. 1995; Jones 1995; Markevitch et al. 1996; Holzapfel et 
al. 1997; Furuzawa et al. 1998). The original Sunyaev-Zel'dovich formula has been derived from a kinetic 
equation for the photon distribution function taking into account the Compton scattering by electrons: the 
Kompaneets equation (Kompaneets 1957; Weymann 1965). The original Kompaneets equation has been 
derived with a nonrelativistic approximation for the electron. However, recent X-ray observations have 
revealed the existence of many high-temperature galaxy clusters (David et al. 1993; Arnaud et al. 1994; 
Markevitch et al. 1994; Markevitch et al. 1996; Holzapfel et al. 1997; Mushotzky & Scharf 1997; Markevitch 
1998). In particular, Tucker et al. (1998) reported the discovery of a galaxy cluster with the electron 
temperature fcsT e = 17.4 ± 2.5 keV. Rephaeli and his collaborator (Rephaeli 1995; Rcphaeli & Yankovitch 
1997) have emphasized the need to take into account the relativistic corrections to the Sunyaev-Zel'dovich 
effect for clusters of galaxies. 

In recent years remarkable progress has been achieved in the theoretical studies of the relativistic 
corrections to the Sunyaev-Zel'dovich effects for clusters of galaxies. Stebbins (1997) generalized the 
Kompaneets equation. Itoh, Kohyama, & Nozawa (1998) have adopted a relativistically covariant formalism 
to describe the Compton scattering process (Berestetskii, Lifshitz, & Pitacvskii 1982; Buchlcr & Yuch 
1976), thereby obtaining higher-order relativistic corrections to the thermal Sunyaev-Zel'dovich effect in 
the form of the Fokkcr-Planck expansion. In their derivation, the scheme to conserve the photon number at 
every stage of the expansion which has been proposed by Challinor & Lasenby (1998) played an essential 
role. The results of Challinor & Lasenby (1998) are in agreement with those of Itoh, Kohyama, & Nozawa 
(1998). The latter results include higher-order expansions. Itoh, Kohyama, & Nozawa (1998) have also 
calculated the collision integral of the Boltzmann equation numerically and have compared the results with 
those obtained by the Fokker-Planck expansion method. They have confirmed that the Fokker-Planck 
expansion method gives an excellent result for ksT e < 15keV, where T e is the electron temperature. For 
ksTe > 15kcV, however, the Fokker-Planck expansion results show nonnegligible deviations from the results 
obtained by the numerical integration of the collision term of the Boltzmann equation. Here it should be 
pointed out that the generalized Kompaneets equation is equivalent to a single-scattering approximation. 
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Thus for high-temperature clusters (fee 7^ >15 keV) the relativistic corrections may underestimate the 
Sunyaev-Zel'dovich effect at high frequencies. 

Nozawa, Itoh, & Kohyama (1998b) have extended their method to the case where the galaxy cluster 
is moving with a peculiar velocity with respect to CMB. They have thereby obtained the relativistic 
corrections to the kincmatical Sunyaev-Zel'dovich effect. Challinor & Lasenby (f999) have confirmed the 
correctness of the result obtained by Nozawa, Itoh, & Kohyama (1998b). Sazonov & Sunyaev (1998a, 
b) have calculated the kinematical Sunyaev-Zel'dovich effect by a different method. Their results arc in 
agreement with those of Nozawa, Itoh, & Kohyama (1998b). The latter authors have given the results of 
the higher-order expansions. 

Itoh, Nozawa, & Kohyama (2000) have also applied their method to the calculation of the relativistic 
corrections to the polarization Sunyaev-Zel'dovich effect (Sunyaev & Zcl'dovich 1980b, 1981). They 
have thereby confirmed the result of Challinor, Ford, & Lasenby (1999) which has been obtained with a 
completely different method. Recent works on the polarization Sunyaev-Zel'dovich effect include Audit & 
Simons (1998), Hansen & Lilje (1999), and Sazonov & Sunyaev (1999). 

In the present paper we address ourselves to the numerical calculation of the relativistic corrections to 
the thermal Sunyaev-Zel'dovich effect. As stated above, Itoh, Kohyama, & Nozawa (1998) have carried out 
the numerical integration of the collision term of the Boltzmann equation. This method produces the exact 
results without the power series expansion approximation. In view of the recent discovery of an extremely 
high temperature galaxy cluster with ksT e = 17 .4 ± 2.5keV (Tucker et al. 1998), it would be extremely 
useful to present the results of the numerical integration of the collision term of the Boltzmann equation in 
the form of an accurate analytic fitting formula. 

Sazonov & Sunyaev (1998a, b) have reported the results of the Monte Carlo calculations on the 
relativistic corrections to the Sunyaev-Zel'dovich effect. In Sazonov & Sunyaev (1998b), a numerical table 
which summarizes the results of the Monte Carlo calculations has been presented. This table is of great 
value when one wishes to calculate the relativistic corrections to the Sunyaev-Zel'dovich effect for galaxy 
clusters of extremely high temperatures. Accurate analytic fitting formulae would be still more convenient to 
use for the observers who wish to analyze the galaxy clusters with extremely high temperatures. This is the 
motivation of the present paper. For the analyses of the galaxy clusters with extremely high temperatures, 
the results of the calculation of the relativistic thermal bremsstrahlung Gaunt factor (Nozawa, Itoh, & 
Kohyama 1998a) and their accurate analytic fitting formulae (Itoh et al. 2000) will be useful. 
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The present paper is organized as follows. In § 2 we give the method of the calculation. In § 3 we give 
the analytic fitting formula. Concluding remarks will be given in § 4. 



2. BOLTZMANN EQUATION 

We will formulate the kinetic equation for the photon distribution function using a relativistically 
covariant formalism (Berestetskii, Lifshitz, & Pitaevskii 1982; Buchler & Yueh 1976). As a reference system, 
we choose the system which is fixed to the center of mass of the cluster of galaxies. This choice of the 
reference system affords us to carry out all the calculations in the most straightforward way. We will use 
the invariant amplitude for the Compton scattering as given by Berestetskii, Lifshitz, & Pitaevskii (1982) 
and by Buchler & Yueh (1976). 

The time evolution of the photon distribution function n(u) is written as 

= ~ 2 / Wf d3p ' d3k ' w {n{uj)[1 + "WVW n{J)[l + "Mtf^')} > 

2ujuj'EE' ' 1 ' ' 

* - - 7) M;^)" -*»'(; + ;?). ( "> 

k = -2(p- k) = -2uE ^1 - i^icosa^ , (2.4) 

k' = 2(p- k') = 2uj'E ^1 - ^icosa'^) . (2.5) 

In the above W is the transition probability corresponding to the Compton scattering. The four-momenta 
of the initial electron and photon arc p = (E,p) and k = (ui,k), respectively. The four- momenta of the final 
electron and photon are p' = (E' ,p l ) and k! = (u>',k'), respectively. The angles a and a' are the angles 
between p and k, and between p and k', respectively. Throughout this paper, we use the natural unit 
h = c = 1 unit, unless otherwise stated explicitly. 

By ignoring the degeneracy effects, we have the relativistic Maxwellian distribution for electrons with 
temperature T e as follows 



f(E) = 



e {(B-m)-( ( u-m)}/fe B T e + j 
e -{K-(p-m)}/k B T c ^ ^2.6) 



where K = (E — m) is the kinetic energy of the initial electron, and (/i — m) is the non-relativistic chemical 
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potential of the electron. We now introduce the quantities 



k B T e 



(2.7) 



Substituting equations (2.6) - (2.8) into equation (2.1), wc obtain 

^InT = ~ 2 J (2^ rfVrf3fc ' W f(E) { [1 + "(^"M - t 1 + n(w)]n(w')e Al } • (2-9) 

Equation (2.9) is our basic equation. We will denote the Thomson scattering cross section by <tt, and the 
electron number density by N e . We will define 

« e = ^, (2.10) 

Ul e C 



y = a T d£N t 



Jd£N e , (2.11) 



where T e is the electron temperature, and the integral in equation (2.11) is over the path length of the 
galaxy cluster. By introducing the initial photon distribution of the CMB radiation which is assumed to be 
Planckian with temperature To 

"»< A 'I = prri' < 2 ' 12 > 
X s *Sv (2 ' 13 » 

We obtain the function F(9 e ,X) by numerical integration of the collision term of the Boltzmann equation 
(2.9). The accuracy of the numerical integration is about 10 . We confirm that the condition of the 
photon number conservation 

J dXX 2 An(X) = (2.15) 
is satisfied with the accuracy better than 1CP 9 . 



we rewrite equation (2.9) as 



We define the distortion of the spectral intensity as 

(2 , 6) 

e x - 1 n (X) 

- y^—rF(6 e ,X). (2.17) 
e — 1 

The graph of F(6 e ,X) is shown in Figure 1. The graph of Al/y is shown in Figure 2. 
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3. ANALYTIC FITTING FORMULA 

We give an accurate analytic fitting formula for the function F(9 e ,X) in equation (2.14) which has 
been obtained by numerical integration of the collision term of the Boltzmann equation. We will give an 
analytic fitting formula for the ranges 0.02 < 9 e < 0.05, < X < 20, which will be sufficient for the analyses 
of the galaxy clusters. For 8 e < 0.02, the results of Itoh, Kohyama, and Nozawa (1998) give sufficiently 
accurate results (the accuracy is generally better than 1%). 

We express the fitting formula for 0.02 < 6 e < 0.05 as follows: 
An(X) 



n (X) 



yF(0 e ,X) 



6 e X 



X 



= y i prry ( y o + 6 ^ + °i Y ^ + + KY A ) + r \ . (3.1) 



The functions Y , Y\, Y 2 , Y 3 , and Y4 have been obtained by Itoh, Kohyama, and Nozawa (1998) with the 
Fokker-Planck expansion method, and their explicit expressions have been given. 

We define the residual function R in equation (3.1) as follows: 

for < X < 2.5 

n { (3.2) 




where 

6 e = 25(^-0.01), 0.02 < 6» e < 0.05 , (3.3) 
Z = j^(X - 2A) , 2.5<X<20.0. (3.4) 

The coefficients a^j are presented in TABLE 1. The accuracy of the fitting formula for equation (3.1) is 
generally better than 0.1% except for a region of 9 e = 0.05, X > 17, where the error exceeds 1%. 



4. CONCLUDING REMARKS 

We have calculated the relativistic corrections to the thermal Sunyaev-Zel'dovich effect for clusters 
of galaxies by numerical integration of the collision term of the Boltzmann equation. We have presented 
an accurate analytic fitting formula for the thermal Sunyaev-Zel'dovich effect. The fitting formula covers 
all the ranges of the observation of galaxy clusters in the foreseeable future. The accuracy of the fitting 
is generally better than 0.1%. The present results will be useful for the analyses of the galaxy clusters 
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with extremely high temperatures. For galaxy clusters with relatively low temperatures 9 e < 0.02, the 
Fokker-Planck expansion results of Itoh, Kohyama, & Nozawa (1998) will be sufficiently accurate (the 
accuracy is generally better than 1%). 

We thank Professor Y. Oyanagi for allowing us to use the least square fitting program SALS. We 
also thank our anonymous referee for many valuable comments which helped us tremendously in revising 
the manuscript. This work is financially supported in part by the Grant-in- Aid of Japanese Ministry of 
Education, Science, Sports, and Culture under the contract #10640289. 
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Figure Captions 



Fig.l. The function F(6 e ,X). The dotted, dashed, dash-dotted, and solid curves correspond to the 
cases for 9 e = 0.02,0.03,0.04,0.05, respectively. 

Fig. 2. The spectral intensity distortion Al/y as a function of X. The dotted, dashed, dash-dotted, 
and solid curves correspond to the cases for 9 e = 0.02, 0.03, 0.04, 0.05, respectively. 
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Table 1. COEFFICIENTS 



(N 



3=0 



.7=1 



OO 

O - 

■a 

i =14 
O , 
i ^j5 

J .^7 



= 10 



4.212483E- 


3 


-4.444816E-1 


1.073238E+1 


-1.096539E+2 


5.964685E- 


-2 


-1.925063E+3 


3.8769: 


6.658650E- 


2 


3.141377E+0 


-5.060928E+1 


2.931651E+2 


-4.799553E- 


-2 


-1.656944E+3 


9.0227' 


1.064073E4 


-0 


-1.495849E+1 


1.862757E4-2 


-1.026793E+3 


2.276656EJ 


-3 


-1.759309E+3 


1.5129^ 


1.110801E4 


-1 


5.585558E+1 


-3.339812E+2 


1.705455E+3 


-2.884721E- 


-3 


-8.599791E+2 


7.4872 


6.567949E4 


-1 


-2.018768E+2 


1.693840E4-2 


-1.350618E+3 


2.741178EJ 


-3 


-1.865475E+3 


7.0541! 


2.285280E4 


-2 


5.817376E+2 


6.539315E+2 


-9.223397E4-2 


3.300752EJ 


-2 


-3.782007E+3 


5.4132: 


4.882355E4 


-2 


-1.136911E+3 


-1.737435E+3 


2.902696E4-3 


-3.374909E- 


-2 


1.583587E+3 


4.2858: 


6.493386E4 


-2 


1.434664E+3 


1.998541E+3 


-2.732035E+3 


2.384819EH 


-1 


-5.098857E4-3 


2.1306' 


5.247520E4 


-2 


-1.134319E+3 


-1.008644E+3 


1.167943E+2 


3.264981EH 


-3 


-2.006417E+3 


3.9747, 


2.360646E4 


-2 


5.147089E+2 


1.003613E+2 


6.701438E4-2 


-1.143976E+3 


1.427423E+3 


-8.37791 


4.537578E4 


-1 


-1.024535E+2 


4.040647E4-1 


1.039578E+2 


-2.327177E- 


-3 


7.547122E4-3 


-1.0970; 



